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1.0 INTRODUCTION: SCOPE OF WORK AND ORGANIZATION OF REPORT 


The central objective of Phases I, II, and III of this study has been to 
develop mathematically-based computer models which simulate the performance of 
electrically large space-borne imaging radar antennas. Such an approach 
provides the systems designer with a powerful tool for use in the analysis of 
tolerance and environmental effects on the performance of the system, and thus 
greatly reduces developmental and testing costs that would be associated with 
trial-and-error methods of assessing these effects. 

This final report presents the results of the third "^hase of the Antenna Eval- 
uation Study for the Shuttle Imaging Radar (SIR). The principal objectives of 
Phase III were (1) to continue simulations of the effects of excitation errors 
in the performance of a planar array antenna, (2) to apply this model specifi- 
cally, to the SIR-A antenna developed by Ball Brothers, (3) to support measure- 
ments of the SIR-A antenna at JPL by quick-turnaround analyses of the SIR-A 
antenna under a variety of conditions encountered at the JPL West Mesa Antenna 
Range, and (4) to document fully all software developed for NASA Johnson Space 
Center (JSC) during Phases I, II, and III of this contract. In addition, sup- 
port was provided to JSC in the following areas: 


1. Presentation at JSC of significant results concerning other possible 
antenna configurations, potential problems with SIR-A, potential 
problems at X-band, etc. (January 2-3, 1978). A demonstration of 
all software developed for the Tektronix 4051 graphics terminal was 
made at the same time. 

2. Participation in Preliminary Design Reviews (PDRs) and Critical 
Design Reviews (CDRs) of the SIR-A antenna and interface. 

3. Investigation of the effects of reflection of radiated energy from 
the SIR-A antenna in a building similar to the KSC O&C building. 

Some of the tasks originally scheduled for the Phase III effort were either 
performed in Phases I and II (such as the computer model development) or were 
shared with the companion JSC contract for Analysis and Support of SIR-A 
development and test. 


Because of the large amounts of data that can be generated by the computer 
programs developed in all three phases, it seems more logical and efficient to 
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present samples of the data that illustrate the potential and versatility of 
the software rather than provide an ency :lopedic report. Ample references are 
made to other reports so that the interested reader may study the actual data 
generated and conclusions drawn. This final report is divided into the 
following sections: 


Section 2.0 The Shuttle Imaging Radar .Antenna General Footprint Simu- 
lation Program 

Section 3.0 The Shuttle Imaging Radar Antenna Specialized Far-Field 
Pattern Program 

Section 4.0 Ray Tracing Program for the SIR-A Antenna and Payload Bay 
Inside an Enclosed Building 


Section 5.0 Summary 

The examples used to illustrate the use of the computer software are actual 
cases simulated for JSC use. Complete listings of each program may be found 
in the Appendix. All TEKTRONIX software developed in the course of this con- 
tract is available on DC-300A cassette tape compatible with the TEKTRONIX 4051. 
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2.0 THE SHUTTLE IMAGING RADAR ANTENNA GENERAL FOOTPRINT SIMULATION PROGRAM 

The mathematical model and the computer program for antenna simulations were 
developed and tested during Phase I; the majority of the simulation scenarios 
were carried out during Phase II and reported in the Phase II Final Report [2» 
pp. 16-60] and elsewhere [3> 4, 5]. Traditional analyses of the antenna 
footprint on the earth's surface required assumptions of uniform antenna 
illumination, planar geometry, etc., that were too restrictive for the infor- 
mation needed in this study. Hence a more general computer model was devel- 
oped by PSL to gain insight into the effect of antenna surface flatness on: 

1. Beam pointing error 

2. Gain degradation 

3. Main beam spreading and breakup 

4. Sidelobe level degradation 

No changes in the mathematics of the model have been made since Phase I [1, pp. 18-30], 
but some modifications to the computer code have been implemented since Phase 
II [1, pp. 31-37]. These changes will be summarized in this section, 'ind 
detailed information regarding data input and output will be documented. One 
sample program output will be included for completeness and continuity. 

2.1 Modifications to the Footprint Program 

The only modifications to the program which are visible to the user are in the 
auxiliary I/O routines of Figure 1. These routines (Tapel, Tape2, DISK, and 
STATUS) were not implemented at PSL because there was no subsequent demonstra- 
tion of need. These routines should be deleted from Figure 1 (and Figure 8 of 
Reference 1). However, if the need arises, these routines may be added. The 
only potential problem is that of compatibility between computer tapes and 
disk drives. Therefore, before these I/O routines are written, users should 
agree on media types, formats, etc. 
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2.2 User*s Guide to the Footprint Program 

2.2.1 Program Description 

The function of each routine listed in Figure 1 is described. 

MAINPGM “ Performs mostly housekeeping chores and calls other routines 

according to user input commands. It calculates the ’’YAW-TILT- 
TWIST* array, predicted beam center location, plot normalization 
factor (pattern value at predicted beam center location) , 
increments beam pointing angle to cover en;.ire ’’footprint** 
region and initiates profiles along lines of constant latitude 
and longitude intersecting the predicted beam center. 

ANTENNA - Inputs and calculates appropriate antenna parameters: number 

of elements, spacing, phase shift, number of subarrays, polar- 
ization, and element type. After calculations have been per- 
formed, a summary is printed. 

MECH - Calculates mechanical deformation data based on the inputs from 
STRESS, THERML, and MISC. For each subarray, ^ECH calculates 
the average displacement ZAVG, the tilt angles, ALPHAX and 
ALPHAY, and the error coefficient of the bilinear approxima- 
tions. A summary is printed after execution. 

MISC - Inputs any other deforamcion data not covered by STRESS and 
THERML. 

ELEC - Inputs the amplitude, phase, and linear phase gradients of each 
subarray. This information is passed to PAT via conunon block 
ELCTRC, and a summary is printed. 

ORBIT - Inputs orbital parameters and antenna YAW, TILT, and TWIST, 
then prints a summary. 
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PAT * Is the basic pattern subprogram. Given latitude-longltude 

coordinates on the earth’s surface, it performs all translations 
and rotations to determine the beam pointing angle (u,v). It 
then calls AF to find the subarray factor and finally computes 
the sum of all subarrays. 

AF - This function calculates the array factors of a subpanel by 
’’table look-up'’ with linear interpolation if juj < UMAX and 
jvj < VMAX. Otherwise, the AF is computed in the usual way. 

GENER - This routine generates the table used by AF. Table length is 
1001 words: UMAX and VMAX are chosen to include the main beam 
plus the first three side lobes. 

OUTPUT - Translates user commands for printer and pattern output into 
logical switches for the program. 

PROFIL - Is the printer profile plot routine. It will print a one-dimen- 
sional plot down the page for up to 501 data points. The exact 
values of both the abscissa and ordinate are printed with each 
data point, the ordinate also being printed in dB. 

PATCON - Is chiefly a ’’bookkeeping” subroutine for CONTUR. It calls 

subroutine CONTUR three times to generate a complete 151 x 151 
two-dimensional contour map. 

CONTUR - Prints a 51 x 151 contour plot of the footprint. To obtain a 
complete 151 x 151 plot, this routine is called three times by 
subroutine PATCON. Then the three separate plots are pasted 
together to make the composite. 

PLOTl - Is the Calcomp one-dimensional plot routine. It is called 

twice by the main program, and it is used to generate the pro- 
file plots of the principle planes through the main beam. 
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PL0T2 * Is the Calcorop contour map routine. From the data used to gen* 
erate the printer contour map, it draws a continuous contour 
plot of the antenna ground footprint over the rectangular 
region specified by the user. 

PLOTS - Is the Calcomp three-dimensional plot routine. It plots the 
magnitude (in dB) of the footprint pattern over rectangular 
region specified by the user. 

2.2.2 Input/Output Commands and Data Cards 

The program input data has been structured so that user knowledge 
of computer fundamentals is not required. Consequently, the 
program can be used by anyone desiring information concerning 
the footprint of a particular antenna configuration. 

The input has been divided into six modules. These are: 

A. Simulation Information 

B. Antenna Position and Orientation 

C. Antenna Configuration 

D. Antenna Mechanical Parameters 

E. Antenna Electrical Parameters 

F. Output Commands and Parameters 

By including input parameters in these six categories, it is 
possible for the user to simulate a wide range of electrical, 
mechanical, and physical scenarios. Furthermore, a change in 
the simulation does not require that the complete input deck be 
repunched. When related scenarios are being studied, only one 
of the six modules need be changed. 

Simulation Information 


Card 1: Simulation number (15) 

Date (5A4) 
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Card 2: Narative description of simulation (80A1) 

End narrative with ampersand (&) 

Antenna Position and Orientation 


Card 1: Frequency in GHz (FIO) 

Card 2: Orbit altitude (km) (FIO) 

Card 3: Yaw, tilt, twist angles in degrees (3F10) 

Antenna Configuration (Subroutine ANTENA) 

Card 1: Number of subarray sections in aximutb and elevation 

(215) 

Card 2: Number of elements in izirauth and elevation (215) 

Card 3: Element spaceing (cm) in azimuth and elevation (2F10) 

Card 4: Interelement phase shift (degrees) in azimuth and 

elevation (2F10) 

Card 5: Polarization (O=isotropic, l=horizontal, 2=vertical) 

(15) 

Antenna Mechanical Information (Subroutines MECH and MISC) 

Cards 1 Warp array in cm (8F10). The displacement of each 
corner of the subarrays is entered. For example, 
suppose a systematic tilt as shown below is to be 
simulated. Then the data entered would be: 


0.0 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

1.0 

(Card 1) 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

2.0 

3.0 

(Card 2) 

4.0 

5.0 

6.0 

7.0 

8.0 




(Card 3) 


0.0 

1.0 

2.0 

3.0 

4.0 

5.0 

6.0 



1.0 

2.0 

3.0 

4.0 

5.0 

6.0 

7.0 



2.0 

3.0 

4.0 

5.0 

6.0 

7.0 

8.0 
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Antenna Electrical Parameters (Subroutine ELEC) 

Card 1 One card for each subarray section to include: 
Phase shift in azimuth (degrees) 

Phase shift in elevation (degrees) 

Magnitude of excitation (volts) 

Phase of excitation (degrees) 



Output Commands and Parameters 


(Subroutine OUTPUT) 


Card 1: Printer output commands (211) 

PR1=1; Printer profile plots 

PR2=1: Printer contour map of footprint 

Card 2: Plotter output commands (311) 

PL1=1: Plotter profile plots 

PL2=1: Plotter contour map of footprint 

PL3=1: Plotter three-dimensional relief map of 

footprint 

[Cards 3-8: Printer/Plotter output parameters) 

Card 3: Limits on X-axis (degrees longitude) (2F10.0) 

Card 4: Limits on Y-axis (degrees latitude) (2F10.0) 

Card 5: Number of contours to be plotted 

PL2=1: Use Card 5 

PL2=0: Omit Card 5 

Card 6: Contour levels to be plotted (dB) (8F10.0) 

PL2=1: Use Card 6 

PL2=0: Omit Card 6 

Card 7: Limits on printer contour map (dB) (2F10.0) 

PR2=1: Use Card 7 

PR2=0: Omit Card 7 

Card 8: Number of points along each axis for the printer 

contour map (215), maximum of 151 x 151 points 
PR2+PL2+PL3 >1: Use Card 8 


PR2+PL2+PL3 =0: Omit Card 8 
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2.3 Example of Input/Output with the Footprint Program 

Many simulations were run for the seven-panel SIR-A array antennas. The char- 
acteristics of this antenna are listed below: 

Frequency = 1.275 GHz 

Number of Panels = 7 

Element Spacing = 16.75 x 13.00 cm 

Excitation = 1/8, 1/8, 1/6, 1/6, : 6, 1/8, 1/8 


Suppose it is desired to simulate this antenna in earth orbit (altitude - 
200 km) and a tilt angle of 47° with a parabolic bow deformation of one cm 
at the edges. The data card deck would be composed of: 


simulation 

information 


00001 16DEC78 

EXAMPLE 1: BBRC 7 -PANEL .ARRAY WITH NATITIAL TAPER 

AND PARABOLIC BOW OF ONE CM & 


antenna 

position 


1.275 

200.0 

0.0 


47.0 


0.0 


(0000700001 


antenna 

configuration 

< 

0005600016 
16.75 13.00 

0.0 0.0 






i 00001 





mechanical 


/i.O 0.5 

0. 166 

0.0 0.0 0. 

166 0.5 

1.0 

information 


1.0 0.5 

0. 166 

0.0 0.0 0. 

166 0.5 

1.0 



[NOTE; Some blanks have been 

deleted 

in the 



above two 

lines 

to save space 

. . • use 8F10 



1 format] 






'o.o 

0.0 

0.875 

0.0 




0.0 

0.0 

0.875 

0.0 




0.0 

0.0 

1.1667 

0.0 


excitation 


0.0 

0.0 

1.1667 

0.0 




0.0 

0.0 

1.1667 

0.0 




0,0 

0.0 

0.875 

0.0 




0.0 

0.0 

0.875 

0.0 




rii 







111 







-1.0 


1.0 



output 

< 

0.0 

00005 


4.0 





n.o 

3.0 

-10.0 

-20.0 

-27. 



0.0 

-27.0 





' 0015100151 
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The complete output from this simulatioa is shovm on the following pages, with 
the exception that the printer contour plot pages have been pasted together 
and photO' reduced for clarity. The plotter form output has not been included. 
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3.0 THE SHUTTLE IMAGING RADAR ANTENNA SPECIALIZED FAR-FIELD PATTERN PROGRAM 


The footprint program discussed in Section 2.0 is quite general and could be 
used to generate far*field pattern data. However, the computer time consumed 
by the program is Lather large for the far-field results which would be pro- 
duced. For this reason, a separate program was written to quickly compute 
far-fleld pattern profiles from a nominally planar array antenna which has 
been mechanically warped or with electrical excitation errors. Two p ograms 
were written during the course of this study. The first program * titter 

in FORTRAN IV to be run on the IBM 370 computer system at PSL. W -vscution 

time for this program was minimal (less than one minute per profil ^ ..-rated), 
the support time (punching cards, waiting for output, etc.) was excessive. 

Hence a second program was written in BASIC for the TEKTRONIX 4051 Graphics 
System. While the BASIC program takes about ten minutes pec profile, there is 
no support time involved. 

The value of having such a facility as the TEK 4051 can be illustrated by two 
examples . 


1. In March 1978, there was a disagreement between Ball Brothers and 
JPL as to the nominal side lobe level of the SIR-A antenna with a natural 
taper. Using the TEK 4051 and this second program it was possible to resolve 
the question (in favor of Ball) in less than fifteen minutes. 

2. In October/November , 1978, the SIR-A antenna was measured at JPL. 

The preliminary pattern data taken did not agree with previously calculated 
results. The question arose, "Is the measurement in error, or does the antenna 
not work?” PSL ran several simulations independent of Ball Brothers and JPL, 
and suggested that a pseudo near field probe of the individual antenna panel 
radiation be made. Using this data PSL (along with JPL) verified that the 
problems lay in the antenna as the data measured at JPL and the simulations 
run at PSL agreed to within 0.2dB at the -20 dB level. 

The data needed to run either of the two far-field programs is essentially the 
.same as the data needed for the footprint program. The second proj;;. ’.'i! has 
prompted data entry, so a separate section on data I/O need not be included in 
this report. However, a complete example is given below. 
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3 . 1 Example Using the Footprint Program 

The sample program used here is the same as is used in Section 2.3. . ata 

entries are prompted by the computer, and user input has been underline The 
next three pages are copies of computer output from the program. 

Data entry for the FORTRAN program is essentially the same as for the foot- 
print program, with the exception th'it orbit information (altitude, yaw, tilt, 
twist) is not needed for far-field plots. Many examples using both these 
programs were given in the Phase II Final Report [2, Ch3], at the review 
meeting held at JSC on January 2-3, 1978, and in the S.ARTC conference proceed- 
ings [3]. 



tttPSL DEFORHED RHTEHHA FAR-FIELD PATTERN PROGRAH$$» 

ENTER FREQUENCY IN GNZ; 1.275 

ENTER TOTAL NUMBER OF ELEMENTS <NX,NY>: 56i 16 
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SAMPLE OF 

DATA OUTPUT LISTING 


CM o ■T’ o V « ^ ^ r*. 00 rw or» tf) Ok j fw oj fo lo o CO VD lo fo ® N. rw r) fo cj ro 
floooc^osTrooooosroTf iy,rocvi-^iomcsvo'^rvinr>orw®o>'-*a>-*(S>rwc\jjo*xiO 
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^ 000000Q^0CD0 00CD0®GDC5>0CD0Q0<9G>0CD<S>(SDa^CDQQ®0 
2 fonj«4oaiCor^^din^rocvi«-«CDo^cDrwg9in*^ro<vj^Qcrto:)is.vDu*>^ro<vi^o^ 

^ forororooi<wojrucvit^iC^i^^i<vj<\j*^«ii4»-4^^«^<iF^^^^®©CD®CD®®Q©®cD 
w I I I I i I I f I I I I I i I I I I I I I I f I I I I t f t 1 t t 
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4.0 RAY TRACING FOR THE SIR-A ANTENNA AMD PAYLOAD BAY INSIDE AW ENCLOSED BUILDING 

In April 1978 » it was suggested that an end*to*end system test of che SIR-A 
radar and antenna be conducted prior to launch to assess and, if necessary 
correct, system problems arising from electromagnetic Interference (EMI) of 
the radar with other shuttle experiments. Of necessity, this test would have 
had to be conducted inside a closed building (perhaps the 0 & C building at 
Kennedy Space Center). PSL conducted a preliminary study, using approximate 
methods, to assess the feasibility of such a test. The study results were 
communicated to NASA/ JSC by technical memo on 28 April 1978. 

To determine whether additional high level radiation would be reflected by the 
enclosing structure back into the critical payload area, a ray tracing program 
was written for the Tektronix 4051 interactive graphics system. The following 
assumptions were made: 

• The SIR-A antenna pattern would not be perturbed by the surrounding 
environment. 

• The energy emitted from the antenna traveled in ray-like paths 
(geometrical optics approximation) whose directions did not change 
unless incident on a surface (Snell's Law). 

• Diffraction and other higher-order effects were ignored. 

• Reflection coefficients were assumed unity (no energy absorbed by 
the structures) . 

• Power density in the ray direction varied as the inverse square of 
the distance from the antenna (Far-field approximation). 

• Only two dimensions, corresponding to the pallet cross section were 
simulated. Since the SIR-A beam is narrow in the third dimension 
(depth or x-axis), this approximatir ri is fairly good. 

• No external objects (GSE, overhead crane, personnel, et..) were in- 
cluded in the simulations. 

It is not possible to provide accurate figures for reflected radiation from an 
analysis as approximate as this one is; however, it is possible to state that 
significant amounts of energy will be coupled from the SIR-A antenna to other 


shuttle experiments via reflections from surrounding objects. In other words, 
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the test results would be meaningless because a small amount of reflected 
energy would completely mask out the EMI effects of radiation coming directly 
from the antenna. 

4. 1 Use of the Ray Tracing Program 

The program is written for the TEK 40S1; hence, data entry is prompted by the 
computer. Data may be entered in one (or both) of two ways: 

1. By specifying the number of tape file on which data has been pr '^''>usly 
stored 

2. By entering data directly from the keyboard 

After all data has been entered, it may be stored on a tape file by pressing 
User Definable Key #4. Object shapes are entered by specifying the (x,y)- 
coordinates of the endpoints of the line segments making up the object. 

(Several objects may be entered at once by concatenating their respective tape 
files.) In addition, the line segment number representing the antenna and the 
(x,y)-coordinates of the ray source must be entered. For the SIR-A antenna 
mounted in the pallet, this information is given in TABLE I. This data has 
been stored on tape file #2. 

4.2 An Example using the Ray Tracing Program 

Figures 2 and 3 illustrate the problem to be solved. The SIR-A antenna plus 
pallet are placed inside a 50' x 100' rectangle. Data are entered as described 
in Section 4.1, and the ray tracing begins. To stop the trace, user key is 
pressed and the following results are printed: 

R - total length of the ray (in the same units as data input) 

ANGLE - "take off angle of ray (degrees) from horizontal 
X0,Y0 - (x,y) - coordinates where the ray terminates. 
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To restart the program changing only the take-off angle, press user key #3. 

The ray trace program is not limited to this type of problem only. For example, 
it has been used to study the effects of ground surface profiles on reflected 
signal levels for antenna range work. 
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TABLE 1. Boundary surface computer description for EHI hazard study. 

Describing Coordinates 


Surface 

(*0*^0^ 

(xi.yj) 

Descriptioa 

1 

00.000,00.000 

100.000,000000 

Floor 

2 

00.000,50.000 

100.000,50.000 

Ceiling 

3 

00.000,00.000 

00.000,50.000 

Left Wall 

4 

100.000,00.000 

100.000,50.000 

Right Wall 

5 

46.162,04.781 

53.838,4.781 

Pallet 

6 

53.838,04.781 

56.506,07.747 

Pallet 

7 

56.506,07.747 

56.974,11.650 

Pallet 

8 

56.974,11.650 

57.156,12.041 

Pallet 

9 

57.156,12.041 

57.937,12.041 

Pallet 

10 

57.937,12.041 

57.937,12.301 

Pallet 

11 

57.937,12.301 

56.375,12.236 

Pallet 

12 

56.375,12.236 

54.944,08.658 

Pallet 

13 

54.944,08.658 

52.862.6.446 

Pallet 

14 

52.862,06.446 

47.138,6.446 

Pallet 

15 

47.138,06.446 

45.056,08.558 

Pallet 

16 

45.056,08.658 

43.624,12.236 

Pallet 

17 

43.624,12.236 

42.063,12.301 

Pallet 

18 

42.063,12.301 

42.063,12.041 

Pallet 

19 

42.063,12.041 

42.844,12.041 

Pallet 

20 

42.844,12.041 

43.026,11.650 

Pallet 

21 

43.026,11.650 

43.494,07.747 

Pallet 

22 

43.494,07.747 

46.162,04.781 

Pallet 

23 

56.636,12.691 

51.720,17.276 

SIR- A Antenna 

NOTE: 

Additional surfaces may be 
0 & C Building. 

added in order to provide 

more detail of the 


Antenna source location: 54.178,14.9835. 




Figure 2. Reflecting surfaces of the pallet and SIR-A antenna. Surfaces 1 through 
4 describe the four surrounding vails of the O&C Building. 
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Figure 3. SIR-A Antenna and pallet inside a rectangular building. 
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1035.82489938 ANGLE • 48 X0»Y8: 93.2556994843919.8436992762 


35 




36 





1548.799188S1 ANGLE - 46 K8»Y8: 56.4ie2682889f 12*9819387617 
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5.0 SUMMARY 


This report has outlined the work performed during Phase III of the Antenna 
Evaluation Study for the Shuttle Imaging Radar. Greatest emphasis has been 
given in the report to areas which have not been reported elsevtK^re, namely, 
the documentation of all computer software written during all three phases of 
this effort. All three programs which are described in this volume have app] ‘ 
cability to other projects within NASA and elsewhere; hence the user*oriente^ 
format of the report. Technical information regarding algorithms used, mathe* 
matics, programming, etc., may be found in the Final Reports for Phases I and 
II, or in the literature cited at the end of this report. 

If the SIR-A mission is successful, future synthetic aperture radars will un- 
doubtedly be launched into space either as free-flyers or as Shuttle experi- 
ments. Potential users of SAR data will be requiring multiple frequencies, 
multiple swath widths and incidence angles, and dual polarization. For example, 
it has been shown that soil moisture determination is best carried out at about 
4.5 GHz and for incidence angles in the 7®-12^ range, whereas crop classifi- 
cation is optimally accomplished at about 14 GHz and with a 40^-50^ incident 
angle range. When such an experiment is finally implemented, the software 
that has been developed »inder the present effort should be very helpful in 
systems planning and design. 
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